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Barotropic FRW cosmologies are presented from the standpoint of nonrelativis- 
tic supersymmetry. First, we reduce the barotropic FRW system of differential 
equations to simple harmonic oscillator differential equations. Employing the fac- 
torization procedure, the solutions of the latter equations are divided into the two 
classes of bosonic (nonsingular) and fermionic (singular) cosmological solutions. 
We next introduce a coupling parameter denoted by K between the two classes 
of solutions and obtain barotropic cosmologies with dissipative features acting on 
the scale factors and spatial curvature of the universe. The K-extended FRW 
equations in comoving time are presented in explicit form in the low coupling 
regime. The standard barotropic FRW cosmologies correspond to the dissipa- 
tionless limit K = 0. 

Comoving FRW barotropy 

Barotropic FRW cosmologies in comoving time t obey the Einstein-Friedmann 
dynamical equations for the scale factor a(t) of the universe supplemented 
by the (barotropic) equation of state of the cosmological fluid 



where p and p are the energy density and the pressure, respectively, of the 
perfect fluid of which a classical universe is usually assumed to be made of, 
k = 0,±1 is the curvature index of the flat, closed, open universe, respec- 
tively, and 7 is the constant adiabatic index of the cosmological fluid. 
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Bosonic FRW barotropy 

Passing to the conformal time variable rj, defined through dt = a(?])dry, one 
can combine the three equations in a single Riccati equation for the Hubble 
parameter H (r/) (we shall use either gj- or ' for the derivative with respect 
to r] in the following) 

Hq + cHg + kc = , (1) 

where c = §7 — 1. 

Employing now Ho(?7) = -— one gets the very simple (harmonic oscilla- 
tor) second order differential equation 

w" - c • c Kjb w = , (2) 

where c Ki b = —kc. Moreover, the particular Riccati solutions H^ = — tancr; 
and Hq = cothcr/ for k = ±1, respectively, are related to the common fac- 
torizations of the equation @ 

A + cH ) (A - cH ) w = w" - c(H + cH§)w = . (3) 

Borrowing a terminology from super symmetric quantum mechanics, we 
call the solutions w as bosonic zero modes. They are the following. For 

K = l 

wi,b ~ cos(ct? + d) -> ai >b (r/) ~ wj /c , 
where d is an arbitrary phase, whereas for k = — 1 one gets 

w_i jb ~ sinh(c?7) -> a_i >b (r?) ~ . 



Fermionic FRW barotropy 

A class of barotropic FRW cosmologies with inverse scale factors with re- 
spect to the bosonic ones can be obtained by considering the supersymmetric 
partner (or fermionic) equation of Eq. (J3J) which is obtained by applying the 
factorization brackets in reverse order 

^ - cH ^) (A + C H ^) w = w" - c(-H + cH§)w = . (4) 
Thus, one can write 

w - c • c K m = , (5) 



where 



r (ri \ _ H ' , p tt2 _ J c(l + 2tan 2 c?7) if k = 1 
c Kif (r?)--H + cH - c( _ 1 + 2coth 2 cr?) ■ liK = _ 1 
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denotes the supersymmetric partner adiabatic index of fermionic type as- 
sociated through the mathematical scheme to the constant bosonic index. 
Notice that the fermionic adiabatic index is time dependent. The fermionic 
w solutions are 

w i>f = 1 C , a\ ~^ a i,f(^) ~ [cos(c?? + d)]~ 1/c , 

cos (cry + d) 

and 

w -!. f = ■ U ^ a-yW ~ [sinh(cr7)] _1 / c , 

smn(cr/J 

for k = 1 and k = — 1, respectively. 

We can see that the bosonic and fermionic barotropic cosmologies are 
reciprocal to each other, in the sense that 

a ±,b a ±,f = const . 

Thus, bosonic expansion corresponds to fermionic contraction and viceversa. 

Uncoupled fermionic and bosonic FRW barotropies 

A matrix formulation of the previous results is possible as follows. Introduc- 
ing the following two Pauli matrices 

a = -i<7 y = -i ^ *? ' ^ and (3 = a x = ^ J J^J 

we can write a cosmological matrix equation 

a y B v W + <r x (icH )W = , (6) 

where W = ( Wl | is a two component 'zero-mass' spinor. This is equiva- 

V ™ 2 J 

lent to the following decoupled equations 

D^wi + cH u>i = (7) 
-B v w 2 + cH u> 2 = . (8) 

Solving these equations one gets w\ oc 1/ cos(cr/) and W2 oc cos(cr/) for k = 1 
cosmologies and w\ oc l/sinh(cr?) and W2 oc sinh(cr7) for k = —1 cosmologies. 
Thus, we obtain 

7 Wl \ = ( Wf \ 

W 2 I [ W b I 



w 



This shows that the matrix equation contains the two reciprocal barotropic 
cosmologies on the same footing as the two components of the spinor W. 
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Coupled fermionic and bosonic cosmological barotropies 

There is a simple way to couple the two spinorial components by means of 
a constant parameter K. Indeed, we write 



[a y D„ + CT x (icH + K)]W = KW , 



0) 



where K is equivalent to the mass parameter of a Dirac spinor. Eq. © is 
equivalent to the following system of coupled equations 



D^wi + (icHo + K)wi = Kw2 
-D r? w 2 + (icH + K)w 2 = Kw 1 . 



(10) 
(11) 



These two coupled first-order equations are equivalent to second order dif- 
ferential equations for each of the two spinor components. 
The fermionic spinor component can be found directly as solutions of 



Wi 



c [ c i,i{v) + 2iKtancr7]w^ = 



1 w 1 — cfc^f (??) — 2iKcoth ci]]w 1 =0 
whereas the bosonic components are solutions of 



Wo 



c[c — 2iKtanc??]w2" = 



w 2 + c[— c + 2iKcoth ct?]w 2 = 



for k = 1 
for n = —1 



for k = 1 
for k = — 1 



(12) 



(13) 



The solutions of the bosonic equations are expressed in terms of the Gauss 
hypergeometric functions 2F1 of complex parameters that can be written in 
explicit form : 



Z2 k2 w+(r ? )=^zf 2 F 1 



ki 



ki + k 2 + l,ki +k 2 ,l + 2ki ; 



.Be-Ki+afciVf I)* 2Fl 



and 



z 4 k % 2 (r,)=Cz^ 2 F 1 



k 3 



-ki + k 2 ,-ki+k 2 + l,l-2ki; 



k 3 + k4,k3 + k4 + l,l + 2k 3 ;^ 



Z3 



Zl 

2 
zi 



-k 3 + k 4 + 1, -k 3 + k 4j 1 - 2k 3 ; -j 



(14) 



(15) 



where the variables z, (i = 1, ...,4) are given in the following form: 

zi = itan(c?7) — 1, z 2 = itan(cr/) + 1, z 3 = coth(cr/) + 1, z 4 = coth(cr/) — 1, 
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respectively. The k parameters are the following: 

1 / 2K\5 1 / 2K\i 

^'^i 1 --) ■ k2 = 2( 1 + T 

and 

1 / 2K\5 1 / 2K\5 

k» = 5 (i+.-) , k 4 = -(i-i_ ^ 

whereas A, B, C, V are superposition constants. Plots of these modes are 
given in Ref. 2b. 

Based on these K zero- modes, we can introduce bosonic scale factors and 
Hubble parameters depending on the parameter K 

a K ,+ = (w^ , H+(77) = i(logw+)' (16) 

and 

&K,- = (w 2 ~) 1/c , Hk(t?) = ~ (logwa)' , (17) 

and similarly for the fermionic components by changing w^ to w^ in eqs. (|16|) 
and l|17|). respectively. 

Comoving K-coupled FRW barotropy: Small K regime 

Introducing the notations A K = -K(^) K =o and F K (t) = (l + |^J, one 
can show that in the small K/c limit the comoving time equations can be 
written as follows: 



a K 



( F i(t)-s)-fe=-y^+3p), 

^Sk^ 2 Fi(t) = 8nGp - (K+2 "k } 

p = ( 7 - l)p . 



3 aX 



(ax could be either aK,+ or aK,- depending on the k case we take into 
account). 

Interpretation 

We come now to the interpretational issue. We consider only the small K 
regime as realistic. Then, the effects of K show up only on the geometri- 
cal quantities without any change in the barotropic equation of state. The 
parameter K introduces an imaginary part in the cosmological Hubble pa- 
rameter H. Since the latter is the logarithmic derivative of the scale factor of 
the universe one comes to the conclusion that the supersymmetric techniques 
presented here are a supersymmetric way to take into account dissipation 
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and instabilities of barotropic FRW cosmologies. However, since K is also a 
coupling parameter between fermion and boson components, the dissipation 
and instabilities belong to the cosmological epochs that occurred before the 
supersymmetry breaking. 
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